Mon. Not. R. Astron. Soc. 000, [THIS] (2002) Printed 16 December 2010 (MN 1^1^ style file v2. 2) 



Inclination-angle of the outflow in IRAS 05553+1631: 
A method to correct the projection effect 

Shuo Kong^*, Yuefang Wu^f 

^School of Physics, Peking University, Beijing 100871, China. 

^Department of Astronomy, School of Physics, Peking University, Beijing 100871, China. 
Accepted by the Monthly Notices of the Royal Astronomical Society 



ABSTRACT 

A mapping study of IRAS 05553+1631 was performed with ^^CO J=3-2 and 
^■^CO i—2-\ lines observed by the KOSMA 3 m-telescope. A core with a size 
of 0.65 pc and with a LTE mass of 120 Mq was defined by the mapping 
with i^CO J=2-l Une. We have identified a bipolar outflow with ^^CO J=3-2. 
For accuracy in the calculation of outflow parameters, overcoming the projec- 
tion effect is important. We propose a new method to directly calculate the 
inclination-angle 6. We establish two basic equations with the help of outflow 
contour diagram and finally obtain the " angle function" and the " angle equa- 
tion" to derive 0. We apply our method to the outflow of IRAS 05553+1631, 
finding that Ouue is 73° and Ored is 78°. Compared to the parameters initially 
estimated under an assumption of 45° inclination-angle, the newly derived 
parameters are changed with different factors. For instance, the timescales for 
the blue and the red lobes are reduced by 0.31 and 0.21, respectively. Larger 
influences apply to mechanical luminosity, driving force, and mass-loss rate. 
The comparisons between parameters before and after the correction show 
that the effect of the inclination- angle cannot be neglected. 
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1 INTRODUCTION 



Massive star formation (MSF) has attracted much attention. It has enormous impact on 
the natal interstellar medium (ISM) and on the evolution of Galaxy. MSF is believed to 
originate in dense molecular cores (DMCs) that can be traced by CO and its isotopologues. 
With the help of mapping, we can reveal the structure of molecular cores and investigate 
the MSF activity taking place within them. 



H igh velocity outflows ha ve also been intensively studied. First uncovered in 1976 ( IZuckerman. Kuiper fc 



1976 



Kwan fc Scoville 



19761), mo 



have attracted much attention (IWu et al. 



ecular outflows toward massive young stellar objects (YSOs) 



2004 ). Perhaps tracing the earliest stage (ILada 



19851 ) of star formation, molecular outflows are critical in the debate about two mechanisms 



for MSF: massive stars forming through accretion-disk-out flow or via collision-coalescence 



(IWolflre fc Cassinelli 



1987 



Bonnell. Bate fc Zinnecker 



19981 ). The properties of outflows re- 



veal the mass-loss phase before the main sequence. What we observe in the sky is not the real 
outflow but its two-dimensional projection. Therefore, it is essential to know the inclination- 
angle 6 between the outflow axis and our line-of-sight. For instance, it will introduce a factor 
of 1/ cos^^ to the momentum. So far the determination of the inclination-angle has been some- 
what neglected though s ome previous resear chers have made a few attempts by modeling 



(ICabrit fc Bertout 



1986 : 



Meixner et al. 



20021 ). 



Until now, authors usually took assumpti ons for the inclinat ion-angle. For instance. 



Goldsmith et al. 



(119841 ) hypothesized 45° while 



Garden et al. 



(|l99lh supposed 60°. The out- 



flow parameters derived from these assumptions can only be meaningful statistically, but 



they are not accurate. In addition, wor 



flows (IBally &: Lada 



1983 



i has been done to study the collimation of out- 



Wu fc Huang 



19981 ). However, the projection effect should also 



be included. 



In order to directly calculate the inclination-angle and to eliminate or reduce outflow 
parameter uncertainties, we propose a new method in this paper. In the next section, we 
will describe our observations of IRAS 05553+1631. In section[3l the results will be presented. 
We will introduce our method and discuss some of its properties in section HI Section [5] is a 
brief summary of our work. 
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2 OBSERVATION 

The observations were performed with the KOSMA 3 m submilhmeter telescope at the 
Gornergrat Observatory in Switzerland in 2002 to 2004. The dual-channel-SIS receiver was 
used. We observed simultaneously transitions in the 230 GHz tuning range (^^CO J=2-l) 
and in the 345 GHz tuning range (^^CO J=3-2). The beam sizes were 120" at 230 GHz and 
80" at 345 GHz. The beam efficiencies for 230 GHz is 0.68 and for 345GHz is 0.72. The 
integration time for each point was about 1.5 min for both ^^CO J=3-2 and ^^CO J=2-l. 
The system temperatures were about 160 K for ^^CO J=2-l and 235 K for ^^CO J=3-2. 
The spectral resolutions were 0.22 km/s and 0.29 km/s for ^^CO J=2-l and ^^CO J=3-2, 
respectively. We used the position-switch mode in the observations. The map step size is 1 
arcmin. 



3 RESULTS 

Both ^'^CO J=2-l and ^^CO J=3-2 emission were detected and mapped. The typical ^^CO 
spectrum is shown in the left panel of Figure [H There is a single component in this molec- 
ular region. Table [1] summarizes the parameters. Column (1) is the IRAS source name. 
Columns (2) to (5) are the positions of the source, including both the 1950 coordinates 
and the 2000 coordinates. Columns (6) to (8) list the ^^CO J=2-l spectrum parameters 
derived from Gaussian fitting. Co l umns (9) to (11) exhibit the IRAS flux parameters. 
According to Iwood fc Churchwelll ( 119891 ). IRAS sources that satisfy Log(F25/Fi2)>0.57, 
Log(Feo/Fi2)>1.30 and that peak at 100 /im are promising candidates for UC Hn regions. 
From Table [1] we can see that the flrst two criteria are satisfled. We conflrm that IRAS 
05553-M631 has its maximum flux at 100 fim. The results show that IRAS 05553+1631 is 
a candidate for UC H n region. The right panel of flgure [T] shows the core contours. The 
grey-scale background is MSX 8.28 fim emission. The three MSX sources are named as Si, 
S2, and S3, respectively. Spatially, Si and S2 are both possibly associated with the core. 
We also made use of the Two Micron All Sky Survey (2MASS) data. Si has one 2MASS 
counterpart 05581473-M631070 whose magnitudes in the J, H, and bands are 13.103, 
10.752, and 9.041, respectively; S2 has one 2MASS counterpart 05581574-1-1631373 whose 
magnitudes in the J, H, and bands are 12.571, 11.886, and 11.565, respectively. Si is the 
reddest source and more likely to be correlated with the molecular core. 

Figure [2] presents the ^^CO J=3-2 results. In the P-V diagrams (c) and (d), the wide 
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wings are obvious. The velocity ranges of the high-velocity gas were determined as -8 km 
s"^ to 3.1 km s~^ and 8.5 km s~^ to 13 km s~^ (hereafter refered to as HV-ranges) for blue 
and red lobes, respectively. The HV-ranges were integrated into the outflow contour diagram 
(Figure [3]). The outermost blue and red contours are 50 per ce nt of the blue and red peak 
intensities, respectively. The contours show a bipolar structure. IShepherd et al.l (Il998[ ) has 
identified YSO G192.16 as the driving source. Both their work (see Figure 1 of their paper) 
and our result (Figure [3]) show that the blue lobe is more extended from G192.16 than the 
red one. G192.16 is in good alignment between the two outflow lobes, and we adopt G192.16 
as the driving centre. 



The core parameters are summarized in Ta ble El Column (2) lis ts the 1.4 kpc distance 



calcu^ 



ated using the Galactic rotation model (iBrand fc Blitzl Il993l ) while IShepherd et al 



(Il998l ) gave a distance of 2.0 kpc. Column (3) gives the velocity range for the ^'^CO J=2-l 
line. Column (4) displays the core radius which is half the mean value of two linear lengths 
of the half maximum contour. Columns (5) to (8) are derived from the radiation transfer 
equation. A ratio of 8.9x10^ for [H2/^^C0] was adopted to calculate the H2 column density 
(column 8). In column (9), we give the H2 density calculated from the column density and 
the diameter (which is twice the radius). In the calculation of LTE mass which is show n in 



column (10), the mean atomic weight factor of 1.36 w as introduced (Garden et al.lll99ll ). In 



Estalella et al. 



I993I). 



column (11), the virial mass was calculated following 

Table [3] summarizes the outflow parameters. Column (1) indicates the outflow lobes. 
Column (2) lists the outflow sizes. The thick contour was covered by an ellipse whose major 
axis is the linear length from west to east and whose minor axis is the linear length from 
north to south. The outflow sizes are the lengths of semimajor axes. Column (3) lists the 
masses for blueshifted and redshifted gas, respectively. The beam averaged column density 
for the ^^CO J=3-2 transition is 



1, (Te^ + 0.92) 
4.8 X 10^2 ^ ' 



T 



rdw 



(1) 



exj9(-33.2/Tea;) J r/^ [1 - exp(-r)] 
where Te^ is the core excitation temperature. The integration is ove r the HV-range. The ^^CO 



J=3- 2 transition is usually optically thick even in the lin e wings (IPlanabeck, Snell fc Loren 



1991 



For the mass 



19831 ). therefore we adopted a mean optical depth r of 4 ( iGarden et al. 
calculation, an [H2/CO] ratio of 10^ was taken. Momentum and energy are displayed in 
columns (4) and (5). The adopted velocity v is the line-of-sight velocity divided by cos^ 
where 9 is the inclination-angle. Column (6) shows the dynamical timescales which were 
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estimated as R/f where R is the distance from the driving centre to the peak of high- velocity 
gas. The mechanical luminosity, driving force, and mass-loss rate are listed in co lumns (7) 
to (9 ), respectively. They were calculated as E/t, P/t, and V jitYw), respectively ( iWu et al. 
20051 ). where is the wind velocity assumed to be 100 km/s. In the calculation, we first 
assumed 6'=45°, the calculation will be improved with our method (see section H72|) . 



4 DISCUSSION 



4.1 Core 



Table [2] shows that the core has a size ~0.65 pc and that VIlte is equal to 1 20 Mf^. The 
^^CO line width is 2.5 km s~^ which is larger than those of low-mass cores ( iMyers et al. 



19831 ). These results suggest that the core is massive. According to virial theory, a core is 
gravitationally bound when its mass is larger than its virial mass. In our results, this is 
not satisfied, which is possibly due to the [H2/^'^C0] ratio adopted or the optically thick 
assumption for ^'^CO J=2-l emission; alternatively, the core may be unbound. 



4.2 Outflow 

In our calculation for outflow parameters, 45° was assumed as the inclination-angle. Actually, 
almost all the outflow parameters depend on the angle. For instance, the velocity of outflow 
gas is the projection along our line-of-sight. Therefore a factor of 1/ cos^^ should be introduced. 
The same goes for momentum, timescale, etc., but sometimes with different factors. The 
inclination-angle also affects the coUimation factor. Thus, accurat e determination of 9 is 



important. Unfortunately, until no w, people assume a value for 6 (iGoldsmith et al 



Garden et al. 



1991 



ZhuetaL 



1984 : 



20061 ) . In the following, we propose a new method to improve 



this situation. 



4-2.1 The inclination- angle 9 

Our method to derive 6 requires knowledge of the outflow geometry, so before the calculation, 
it is necessary to discuss something about the shapes of outflows. Generally, without impact 
from the surrounding materials, the molecular outflows are likely to be a xially symmetric 



his c ould be in ferred frorn the s ymmetry of the outflow contour diagrams (IGoldsmith et al. 



(11984), Figure 2: 



Wu et al. 



(l2005h. Figure 5 c: 



from modeling ( ISnell. Loren fc Plambeckl Il980l : 



Ballv fc Lada f 



Meixner et al. 



1983 ), Figure 13 NGC 2071) and 



2002 



Meyers-Rice fc Lada 
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19911 ). In our results, the P-V diagram (Figure |2]) and the outflow contours (Figure [3]) 



show good symmetry, especially the blue lobe. In addition, the ellipse-like outflow contours 
indicate a projection of a conical outflow (see Section r4.2.3p . Since a cone is an ideal form 
for an axially symmetric outflow (see Section |4. 2. 2[ Paragraph 2), we assume that the shape 
of the outflow is a cone. We also assume that high velocity components form cone-shell like 
structure whose emission intensity decreases from inside to outside. The cone could be called 
a "cone onion". 

When we observe, what we see is the outflow projection on the sky. Thus, if the outflow 
axis is at an angle 9 to our line-of-sight, the contours will appear as ellipses (Figure lU 
the contours could also be parabola or hyperbola, see section H:.2.3p . Figure |5] presents the 
calculation of the inclination-angle 9. In Figure [5](a), one can see the driving centre and an 
elliptical outflow contour. From Figure [U^b), a relation of the four angles: (3\, /32, 9i, and 9 
can be expressed by equation ([2]) 

tan 6*1 tan/3i f3i , . 

= — « — (2) 

tan 9 tan /32 /32 

where Pi and P2 are the angle distances in the contour-diagram (/3i accounts for the angle 
between driving centre and point D; P2 is the angle between driving centre and axis point A). 
Pi and P2 are usually very small (about several arcmin), thus the small angle approximation 
is valid. 9i and 9 are unknown. Then, we have to flnd another equation for either 9i or 9. 
Notice that the length of line AB does not change as a function of 9. In Figure [5](a), we have 



AB = AO tan a (3) 

where a is the angle between OA and OB. Overlines indicate length. If the outflow contour is 
not symmetric about line AO, we take a as half the angle ZBOC. We have already assumed 
the shape of cone, so the opening angle 9-9i is the same for the entire cone. Therefore in 
Figure E](b) we have 



AOi = A02sm9 (4) 



AB = A02tan{9 - 9i) (5) 



where Oi and O2 are the same points as O in the plane of the sky. AOi is equal to AO. 



Thus, combining equations (E]), (H]), (E]), and substitute AOi with AO we have 

tan(g - 9i) 

tana = (dj 

sin 9 

Putting together equation ([2]) and equation (Q, we can solve the inclination-angle 9 by 
eliminating 9i. Since P2, and a can be derived from the outflow contour, we try to 
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further manipulate the two equations to obtain a better expression. We define the "angle 
constants" 

P=^i (7) 

and 

Q = tana (8) 

Combining equations ([2]) and (Q and eliminating 6i, we have 

1 - P 1 

— - Ptan^tan^ = 1 (9) 

Q cosO 

Defining the "angle function" 

1 - P 1 

A(e) = — - P tan e tan 6, (10) 

Q cos 9 

after determining P and Q, we obtain the inclination-angle by solving the "angle equation" 

m = 1 (11) 

4-2.2 Comparison 

We utilize our method on IRAS 05553+1631. Figure [6] shows the derivation process. In 
general, the shapes and positions in Figure [6] are similar to those in our model (Figure 
[5]^a)). Table m presents the comparison between newly derived and old outfiow parameters. 
We find the inclination-angles to be 73° and 78° for blue and red lobes, respectively. This 
suggests that the two outfiow lobes are in good alignment. In Table HJ compared with the 
former results, the newly derived parameters have rather large corrections. For the blue 
lobe, the momentum is enlarged by a factor of (cos45°/cos73° ~) 2.4, as the velocity is 
corrected by 1/cos^^. The correction for energy is proportional to that of which is about 
5.8. The timescale is reduced by a factor of (cot73°/cot45° ~) 0.31. The correction factors 
for mechanical luminosity, driving force, and mass-loss rate are 19, 7.7, and 7.7, respectively. 
For the red lobe, the analysis is similar, but the factors are different. These results suggest 
that the correction from inclination-angle cannot be neglected and our method is practical. 

From the process of our method one can see that the cone assumption ensure the viability 
of equation ([6]). If not, the equation would be an approximation. Let 6o=6-6i, the Taylor 
expansion of tan^ around 6o is (keeping the first order) tan6'=tan^o+ (^coseo f • -^^^ hhie 
outfiow lobe in our case, tan6'o=0.143, if A^^=l° (0.017 radians), then the second term of 
the Taylor expansion would be ~0.018, the uncertainty is about 12 per cent. For the red 
outfiow lobe, the uncertainty is about 4 per cent. 
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4-2.3 More properties 



Geometrically, when a cone is cut by a plane, the resulting curve can be an ellipse (including 
a circle), a parabola, or a hyperbola, depending on the angle between the cone-axis and the 
plane. For an ellipse, the larger the inclination-angle is, the larger its eccentricity will be, 
and vice versa. When 6 approaches 90° we obtain a parabola or a hyperbola. Theoretically, 
all of the curves can be utilized for the calculation. But the ellipse has the advantage that 
it is usually much more conspicuous and easy to manipulate. The parabola and the hyper- 
bola, on the contrary, can be confusing because of the large inclination-angle. The outflow 
edge from the driving source has very weak emission and the contours exhibit a fan-like 
structure. The centre and the curve will be ambiguous. The opposite extreme is when the 
inclination-angle becomes zero and the ellipse changes to a circle. Fortunately, the angle 
equation ffTTj) still works. In the angle function f lTOj) . the coefficient of the first term 



Q 



is ( AB^^ ~ '~^\r^ ~^)- ^s ^ approaches zero, the projection shape becomes a circle and 



M AB AB 

OA 



= — Meanwhile, cos^ ^1 and the second term vanishes. Thus K{6 — >0)^1. The angle 



AB 



equation ffTTl) is tenable in this extreme case. However, when 9 approaches zero, use of this 
method could be problematic if the resolution is low. Another property is the coUimation. 
In our view, it can be manifested by the opening angle 9-9i. Larger 9-9i means lower col- 
limation, and vice versa. Thus it is reasonable to define a coUimation factor as cot(^-^i) 
(=(tanasin^)~^). In IRAS 05553+1631, the factors are 7.0 and 1.6 for blue and red lobes, 
respectively. One thing should be mentioned, the ellipses we used for IRAS 05553+1631 are 
90 per cent contours (Figure [6]) as stronger emission could reduce the errors. If we chose the 
50 per cent contours, the factor would change. Additionally, the method highly depends on 
the accuracy with which the driving centre is located. One question is the identification, and 
the second is the spatial resolution. The identification in volves personal iudgemen t, which 
is accompanied by considerable uncertainty in some cases. ICyganowski et al.l (120091 ) showed 
that 6.7 GHz methanol maser and relevant 24 /im emission usually coincide with each other 
and they are good tracers to the driving centre, this may offer much help for the iden- 
tification. The resolution is usually low for single-dish observations. Using high-resolution 
telescopes and interferometers can improve this. 

In the angle function ( ITOi) . there are two parameters P and Q. Since we can always 
choose an ellipse excluding the driving centre, it is reasonable to assume 0<P<1. Figure 
[7] illustrates the value of Q, which shows the projection on x-z plane. Lines AOi and AD' 
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represent the projections of two conditions for the plane of the sky: (1) AOi perpendicular 
to O2E (solid, hereafter named cl); (2) AD' not perpendicular to O2E (dashed, hereafter 
named c2). Condition c2 is denoted by prime. AE equals the AB of Figure |5l^a). Thus, it is 
easy to see that the maximum of Q is == which is just l/cos(^-^i). Since {9-6i) cannot be 
large and an angle of 60° can only bring Q=2, we consider 0<Q<2. When P=0, the second 
term of the angle function ( [T0|) vanishes and ^=arccos(l/Q). In order to have solution, Q 
must be larger than 1. In fact, in this case Q reaches its maximum (see Figure [7]). The result 
is 0< 9 <71°. When P=l, there is no solution. Figure [H] shows the variation of the angle 
function as a function of P and Q. With a fixed Q, the function decreases as P increases. 
With a fixed P, the function decreases when Q increases. Usually within (0°,90°), there is 
one solution. When Q approaches 1 (Figure EJ^d)), in order to have solution, P must be 
closer to 0. When P approaches 1, the same happens to Q. The variation of angle function 
is sensitive to the values of P and Q. 



5 SUMMARY 

We mapped IRAS 05553+1631 with ^^CO J=3-2 and ^^CO J=2-l lines. A core was identified 
from ^'^CO J=2-l observations. It has a size of 0.65 pc and LTE mass of 120 which is 
lower than the virial mass of 850 Mq. ^^CO J=3-2 mapping revealed a bipolar outflow. Its 
parameters were initially estimated under the assumption of a 45° inclination- angle. A new 
method to directly calculate the inclination-angle 9 was proposed, and was utilized for the 
bipolar outflow of IRAS 05553+1631. We found that 9uue is 73° and 9red is 78°. Parameters 
with the new ^s were compared with the former ones. For the blue lobe, the momentum was 
enlarged from 82 M© km s~^ to 200 M© km s~^ by a factor of 2.4 while the timescale was 
reduced from 8.8x10^ yrs to 2.7x 10'* yrs by a factor of 0.31. The enlarging factors for energy, 
mechanical luminosity, driving force, and mass-loss rate are 5.8, 19, 7.7, and 7.7, respectively. 
For the red lobe, the momentum was enlarged from 11 Mq km s^^ to 36 Mq km s~^ by a 
factor of 3.4 while the timescale was reduced from 6.3x10* yrs to 1.3x10* yrs by a factor of 
0.21. The enlarging factors for energy, mechanical luminosity, driving force, and mass-loss 
rate are 12, 55, 16, and 16, respectively. The results show that a selection of parameters 
were influenced by the inclination- angle 9. 
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Table 1. Observed parameters. 



Source R.aQ 
(1950) 
(1) (2) 


DEC. 
(1950) 
(3) 


R.A. 

(2000) 
(4) 


DEC. 
(2000) 
(5) 


km/s 
(6) 


K 

(7) 


km/s 
(8) 


(9) 


Log{^) 
(10) 


FlOO 

Jy 

(11) 


05553+1631 05 55 18.0 


16 31 00 


05 58 11.5 


16 31 14 


5.5 


3.1 


2.5 


1.70 


2.52 


528 



Columns (2) to (5) present the position. The coordinate years are in parentheses. 
^ System velocity. 

^ Antenna temperature of ^^CO J=2-l spectrum. 
<^ Pull Width at Half Maximum (FWHM). 

^ Columns (9) to (11) list parameters related to IRAS flux density. 



Table 2. Core parameters. 



Source 

(1) 


kpc 
(2) 


RangfQ 

km s~ ^ 
(3) 


RadiuQ 

pc 
(4) 


T,J3 

K 
(5) 


(6) 


JVi£] 

lO^^cm"^ 
(7) 


lO^^cm^^ 
(8) 


n«J3 

10^ cm~^ 
(9) 


(10) 


(11) 


05553+1631 


1.4 


(3.5,7.3) 


0.65 


25.8 


0.22 


6.8 


6.0 


3.0 


120 


850 



°' Distance. 

^ Integrated velocity range in ^^CO J=2-l line. 

Half the mean value of two linear lengths of the thick contour. 

Excitation temperature. 

Optical depth of ^^CO J=2-l line. 
/ Column density of ^^CO. 
9 Column density of H2. 

Density of H2. 
* Mass estimated under LTE assumption. 
i Virial mass l|Estalella et al.lll993h . 
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Table 3. Outflow parameters. 



Lobe 

(1) 


pc 
(2) 


a43 

Mq 

(3) 


Mq km s"-'- 
(4) 


iQ 

10'^'^ ergs 
(5) 


G 

10* yrs 
(6) 


10-1 Lq 
(7) 


10'"- Mq km yr~l 
(8) 


(9) 


blue 


0.65 


6.2 


82 


11 


8.8 


10 


9.3 


9.3 


red 


0.60 


1.9 


11 


0.60 


6.3 


0.80 


1.7 


1.7 



Half the linear horizontal length of the 50 per cent outflow contour. 
^ Outflow mass in the unit of Mq. 
^ Momentum. 

Energy. 

^ Dynamical timescale. The unit is year. 
f Mechanical luminosity. 
9 Driving force. 
Mass-loss rate. 



Table 4. Comparison between newly derived parameters and previous results for IRAS 05553+1631. 
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Figure 1. ^^CO J=2-l spectrum and contour diagram. Tlie contours range from 20 per cent to 90 per cent of tlie peak 
integration in steps of 10 per cent. Tlie thick line shows the half intensity contour. The grey-scale background is the MSX 8.28 
emission. The three MSX sources are denoted Si, S2, and S3, respectively. 
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Figure 2. ^^CO J=3-2 results for IRAS 05553+1631. (a) spectrum, (b) ^^CO J=3-2 grid map. (c)(d) ^^CO P-V diagrams. 
They are along two different longitudes indicated beside the position axes. The dashed lines indicate the start of the blue and 
red lobes, respectively. 
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Figure 3. Outflow contour diagram. The east blue contours represent the blue outflow lobe while the west red contours 
represent the red lobe. The peaks for the blue and red lobe are 14.6 K km s~^ and 4.3 K km s~^. The triangles show the 
positions of the three MSX sources. The square stands for the position of the ^^CO core centre. The black solid circle presents 
the position of G192.16. 



case 1 case 2 




Figure 4. Illustration of the geometry. The cones represent outflows with different cases. Observer is looking up from the 
bottom of the figure. The circles and ellipses are what the observer sees in the plane of the sky in the two cases. Case 1 
represents the outflow axis aligned with the line-of-sight. In case 2, the outflow axis is inclined. It is obvious that in case 2, the 
ellipses are not co-centred. And the projection point of axis is not at the centre of the ellipse. 
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red outflow cone : 
axis and one shell 




Figure 5. Illustrations for the calculation, (a) Outflow contours on the plane of the sky. The blue ellipse on the right represents 
the contour of the blucshifted gas while the red ellipse on the left represents the contour of the redshifted gas. Point O shows 
the position of the driving centre. The xyz-coordinate system is centred on the driving centre, (b) The projection plane 
perpendicular to the y-axis. d is the inclination-angle. Points A and D in (b) are the same as those in (a). Q represents axis 
pointing outward perpendicular to the paper while represents axis pointing inward. 
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Figure 6. Calculation of the inclination-angle. The blue contours on the left represent the blue lobe while the red ones on 
the right represent the red lobe. The black solid circle is the driving centre G192.16. The relevant angles are presented in the 
figure. 



X' 




Figure 7. The illustration of parameter Q. The utilization of letters follows FigureO Condition c2 is denoted by prime. 
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Figure 8. Angle function curves with different P and Q. Each of (a)(b)(c)(d) shows the variation of the angle function with 
the changing of P (from 0.2 to 0.8 with an interval of 0.2) and with a fixed Q ( (a) Q=0.2; (b) Q=0.4; (c) Q=0.6; (d) Q=0.8). 
The 6 axis is in angle system. The blue straight line shows A(9)=l. The d of the intersection point indicates the solution to 
the angle equation. 



